Effect of radiation on transport in graphene 
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We study transport properties of graphene-based p-n junctions irradiated by an electromagnetic 
field (EF). The resonant interaction of propagating quasiparticles with external monochromatic 
radiation opens dynamical gaps in their spectrum, resulting in a strong modification of current- 
voltage characteristics of the junctions. The values of the gaps are proportional to the amplitude 
of EF. We find that the transmission of the quasiparticles in the junctions is determined by the 
tunneling through the gaps, and can be fully suppressed when applying a sufficiently large radiation 
power. However, EF can also generate current, but not only suppress it. We demonstrate that if the 
height of the potential barrier exceeds a half of the photon energy, the directed current (photocurrent) 
flows through the junction without any dc bias voltage applied. Such a photocurrent arises as a 
result of inelastic quasiparticle tunneling assisted by one- or two-photon absorption. We calculate 
current-voltage characteristics of diverse graphene-based junctions and estimate their parameters 
necessary for the experimental observation of the photocurrent and transmission suppression. 
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I. INTRODUCTION 

Since its first fabricatiorP, graphene, a monolayer of 
carbon atoms, has attracted a lot of attention as a can- 
didate for the base material in nanoelectronics. The 
unique feature of this two-dimensional (2D) semiconduc- 
tor is the absence of Jhe gap between the conduction 
and the valence bands^, which allows one to change the 
type of the carriers and vary their density in a wide 
range by applying an external gate voltage. Other out- 
standing properties of the material are the high car- 
rier mobilit y and a submi cron mean free path at room 
temperatur d^ ' ^ ! '^ ' ^ ' ^ ' '^ ' ^ ' ^ ' ^^ '. As a result, a lot of experi- 
mental activity has been devoted recently to investiga- 
tion of possibilities of usin g graphene in the field-effect- 
transistor type application j^ ' ^^ ' ^^'^^ ' ^'^ ' ^^ ' ^^ ' ^^ i 

In order to construct such devices one needs to be able 
to control and switch off currents applying gate poten- 
tials. In conventional semiconductors currents can be 
blocked by applying gate voltages that create sufficiently 
high potential barriers confining electrons in a certain 
region. However, the dynamics of low-energy quasipar- 
ticles in graphen e is d escribed by the zero-mass Dirac- 
type HamiltoniarP^nH^ making possible their reflection- 
less transmission through potential barriers of arbitrary 
strength (the so-called Klein paradox) Such an un- 
usual phenomenon reduces the possibility to switch off 
currents by means of gate electrodes and is a very im- 
portant property of graphene. 

It is also well known for many years that the transport 
properties of diverse semiconducting nanostructures are 
sensitive to an externally applied electromagnetic field 
(EF). Though optoelectronic devices such as radiation- 
controlled field-effect transistors, photodiodes, and light- 
emitting diode s have been f abricated on the basis of car- 
bon nanotubed^^^^^l^, a little attention has been 
paid so far to the possibility to control and switch off 



currents in graphene using external electromagnetic ra- 
diation. 

In this paper we study theoretically the transport in 
graphene subject to a coordinate-dependent potential 
U{t) and irradiated by electromagnetic field. For sim- 
plicity we assume that EF is monochromatic and linearly 
polarized but a proper generalization can be easily made. 
We show that the resonant interaction of propagating 
quasiparticles in graphene with EF leads to the forma- 
tion of a dynamical gap A in the quasiparticle spectrum. 
The value of the gap depends on the intensity S and the 
frequency Q of external radiation. The formation of such 
a dynamical gap is a generic quantum property of sys- 
tems described by the two-bands Hamiltonian, and the 
quantity A/h has the same meaning as the famous Rabi 
frequency for microwave- induced quantum coherent os- 
cillations between two energy levels^^. For conventional 
semiconductors the dynamical gap has been predicted in 
the Ref. |27| Charged impurities in a sample^^, electron- 
phononP^, and electron-electron interactions^^ smear the 
quasiparticle spectrum and wash out the dynamical gap. 
Therefore, the dynamical gap in a uniform semiconduc- 
tor can be observed only if the radiation intensity exceeds 
some critical value Sc Under such conditions it has been 
observed in the optical experiments^^ 

So far only bulk properties (s uch as conduct ivit}!^ 
and coefficient of light absorptioiP^^^^l^ in absence of 
coordinate-dependent potentials have been studied in 
semiconductors with dynamical gaps. In this case an 
experimental observation of a dynamical gap is quite dif- 
ficult and demands considerable efforts. At the same 
time, the question about what happens if the system is 
subject to an external potential has not been addressed 
yet. 

In this paper we consider irradiated graphene junctions 
subject to a non-uniform step-like external potentials, 
such that the energies of conducting electrons are lo- 
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cated inside the dynamical gaps only in some narrow res- 
onant regions. This means that during the transmission 
through the junction each electron in a large energy inter- 
val of the order of potential height has to pass through a 
small "classically forbidden" spatial region, where it has 
to move inside the gap. Thus, the current- volt age char- 
acteristics of the junctions are mainly determined by the 
tunneling through the gaps. 

Changing the values of the gaps (i.e. the intensity 
and the frequency of external radiation) and the heights 
of potential barriers, it is possible to vary the current- 
voltage characteristics, e.g. to suppress Klein tunneling 
in graphene p-n junctions by a sufficiently strong radi- 
ation. A short account of this phenomenon has been 
considered recently in Ref. 1361 

However, it turns out that even more interesting effect 
is possible in such a p-n junction. Here we demonstrate 
that in this system, provided the height of the potential 
is larger than a half of the photon energy, the directed 
current {photocurrent) flows through the junction with- 
out any dc bias voltage applied. Such a photocurrent 
arises as a result of inelastic quasiparticle transmission 
through the junction assisted by one- or two-photon ab- 
sorption. We show, that in the presence of impurities 
and electron-electron interaction in the limit of small ra- 
diation intensities the coefficients of reflection or trans- 
mission at the resonant regions slightly differ from those 
in ballistic case even if the radiation intensity S is much 
below the critical value 5'c, at which the dynamical gap 
vanishes. As a result, for experimentally relevant param- 
eters of the junctions the photocurrent is not affected by 
elastic impurities and electron-electron interactions and, 
we hope, it can be observed rather easily at moderate 
radiation intensities. 

The photocurrent we calculated for a typical graphene 
p-n junctions is by several orders of magni tude la rger 
than those measured in carbon nanotube d^^ * ^^ * ^^ * for 
the same radiation intensities, since the effectively two- 
dimensional graphene junction has accordingly larger 
number of conducting channels. As there is no gap be- 
tween the conduction and the valence bands in graphene, 
and the level of doping can be varied in a wide range us- 
ing gate electrodes, graphene-based photodiode can be 
operated in a wide frequency range of the external radi- 
ation in the far-infrared region. 

The paper is organized as follows. In Secjllj we derive 
a general expression for the current flowing through the 
junction in presence of a monochromatic linearly polar- 
ized radiation. The current is determined by the sym- 
metries of the Hamiltonian and the scattering matrix 
for the system. In order to find explicitly the scatter- 
ing matrix in a spatially inhomogeneous potential in the 
presence of radiation one needs to know the probability 
of tunneling through the dynamical gap, which is calcu- 
lated in Sec. IIIIl Then, 



calculate in Sec. |V| current- voltage characteristics of ir- 
radiated ballistic graphene junctions. In Sec. IVl| using 



the kinetic equation approach, we analyze effects of disor- 
der and electron-electron interaction on the photocurrent 



in graphene p-n junctions. In Sec. |VII| we estimate the 
magnitudes of the photocurrent and the tunneling proba- 
bilities for experimentally achievable radiation intensities 
and parameters of the junctions. 



II. CURRENT- VOLTAGE CHARACTERISTICS 
OF IRRADIATED BALLISTIC GRAPHENE 
JUNCTIONS: GENERIC ANALYSIS 

In this section we derive a general expression for 
the cur rent- volt age characteristics of an irradiated two- 
dimensional graphene strip connected to two reservoirs, 
from now on called the left (L) and the right (R) 
leads. We assume that electrons in the strip move in a 
coordinate-dependent potential U{z)^ varying only along 
the strip, and interact with a time-periodic EF. In the 
present and two subsequent sections we consider a clean 
system, such that the transport in the strip is purely bal- 
listic. 

Another important assumption we use is that far in 
the leads the radiation is either absent or its effect on the 
electron motion and distribution functions is negligible. 
In the next section we show that the latter condition is 
satisfied since the majority of the conducting electrons 
in the reservoirs have their energies not too close to the 
positions of the radiation-induced gaps in their spectra. 

Our derivation is similar to the one presented in the 
Ref. 37 but is adopted for the case of two-dimensional 
graphene strip and leads. Assume, that far in the a-th 
lead the electron states are plane waves, incoming to the 
strip, \e^6y^, or outgoing from it, le,^)^^^ characterized 
only by the energy e and the angle of incidence or 
scattering correspondingly. Note here, that there is no 
intervalley and spin scattering in a ballistic sample and 
thus we disregard these degrees of freedom. 

According to the Floquet theorem^^, the general solu- 
tion of the Schrodinger equation for an electron, moving 
in a static potential and subject to a periodic perturba- 



tion, takes the form ^(t) 



*^t(^), where <^>t(^) is a 



^ ^ in Sec. |rV] we consider electron 

ballistic trajectories in a p-n junction to determine how 
many times and where electrons tunnel during the trans- 
mission through the junction. Using these results, we 



periodic function having the same period T as the per- 
turbation. Therefore, an electron incident with energy s 
on the strip can scatter only into the states with energies 
e^nhn, where Q = 27r/T and n = 0, ±1, ±2, ±3, . . ., 
or, in other words, the particle can gain or lose only an 
integer number of quanta hft. 

Thus, the state |£, scatters into 

I d(l)tai3{e,0;e^nhn,(l))\s^nhn,(l)y^''\ (2.1) 

/3,n 

where we have introduced the amplitude taf3 of scattering 
from the state of an electron incident on the strip from 
the a-th lead at the energy e and angle into the state 
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outgoing from the strip in the /3-th lead at the energy e + 

nhO. and angle (j). The quantity \ta(3{s^ 0;e-\- nhVt^ 

gives the probability to scatter into the angle interval 

Accordingly, 



E 

/3,n 



|tc,/3(£:,6>;£: + nfiQ,0)P = 1. 



(2.2) 



In order to obtain the current flowing through the 
strip we introduce the distribution functions f^'^'^*{e^O) 
of electrons in the incoming/outgoing states of the a-th 
lead. In the incoming states they simply coincide with 
the usual Fermi distribution functions for a given tem- 
perature and voltages applied, and, therefore, they de- 
pend only on the energies £ but not on the angles and 
(j). Using the particle conservation law the distribution 
functions of the outgoing states can be rewritten as 

= IZ / d<l^\tf3a{£^rihn,(l);e,0)\^fj'{e^nhn). (2.3) 

/3,n 

The total current flowing through the strip takes the 
form 

I = ^wJ ^0^evcose{fr{e(p))-fr{e(p),e)XU) 

where W is the junction width, and v is the velocity of 
electrons. The coefficient 4 in the last equation accounts 
for the spin and valley degeneracies of th e quasiparticle 
spectrum in graphene. Using Eqs. (2.2) and (2.3), we 
rewrite Eq. (2.4) as 



,jjr ,'pdpdO 

I — rre^'COS 



(2.n)---^g. 

X {\tL^{e,0',e^nh^l,^)\^fT{e)- 
-\t^L{e + nh^, (/>; e, e)\^ff{e + nh^)) . (2.5) 



Changing the variables in Eq. \2.b\ and using the trans- 
verse momentum conservation law, p*"^(£:*^) sin 6>*^ = 
poutf^^out^ ^[^Qout^ ^^^^ terms with /3 = L 

vanish. Indeed, they are responsible for the backscat- 
tering and therefore cannot contribute to the current 
through the junction. 

If the time-periodic EF possesses the time-reversal 
symmetry {t — t), with respect to some moment of 
time t = 0, e.g. if the EF is sinusoidal, then the relation 

\tLR{e,0,e^nhQ,^)\^ = \tRL{e ^ nhn,^;e,0)\^ (2.6) 

holds, and Eq. ( |2.5[ ) can be reduced to 

' pdpdO 



I = 4W 



I 



rev COS X 



where the function 

PLR{£:e^nhn,0) = J d^\tLR{e,0,e^nhn,^)\^l2.S) 

is the probability for a particle in the state \e^O)Y^ in the 
left lead to be scattered into the right lead into the state 
with the energy e + nhft. 

This is the most general expression for the current that 
can be derived using only the symmetries of the Hamil- 
tonian of the system and the particle conservation law. 
The total current / contains the elastic component (the 
term with n = 0) and inelastic components correspond- 
ing to the terms with n 7^ 0. Let us emphasize that 
the radiation not only induces the current due to the 
inelastic processes but also leads to the modiflcation of 
the probability PLi?(^,^,^) of the elastic scattering. In 
the next sections we show that this probability can be 
strongly suppressed by the external radiation, which re- 
sults in the vanishing of the current through the graphene 
strip. Moreover, in systems where the inelastic scattering 
violates symmetries in such a way that 

PLuie, e + nhQ, 0) ^ Plr^s + nh^, e, 0), (2.9) 

the external radiation generates a photocurrent, i.e. the 
directed current without any external voltage applied to 
the junction^''. 

In a related paper, Ref. ^0 the possibility of the pho- 
tocurrent has been overlooked since there was not consid- 
ered inelastic transmission of electrons through the junc- 



PLR{e, e + nhn, 9) (/f (e) - /^(s + nhn)) ,(2.7) 



tion, corresponding to the terms with n 7^ in Eq. (2.7). 

The probabilities Plr{s, £ + nh^}, 0) can be calculated 
explicitly as soon as the Hamiltonian of the system un- 
der consideration is given. In the next two sections we 
flnd the probabilities Plr{^^ ^ + nWt^ 0) for the quantum 
mechanical problem of electron motion in graphene in an 
external potential U (r) and irradiated by a time-periodic 
EF. 



III. EF-INDUCED DYNAMICAL GAPS AND 
ELECTRON DYNAMICAL TUNNELING 

To elucidate the main phenomena related to the in- 
fluence of radiation on the transport in graphene, we 
flrst consider the modiflcation of the quasiparticle spectra 
due to the presence of a monochromatic electromagnetic 
wave, and their tunneling through smooth potential bar- 
riers. For simplicity, we neglect in the next two sections 
scattering due to impurities and electron-phonon inter- 
action. 

The states of electrons in graphene are conveniently 
described by the four-component wavefunctions, deflned 
on two sublattices and two valleys. Since we consider 
the transport in an inflnite clean sample, we may neglect 
intervalley- and spin scattering and study the propaga- 
tion of electrons for different valleys and spin directions 
separately. 
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Electron motion in the time-dependent EF is described 
by the non-stationary Schrodinger equation, 



(3.1) 



where H is the fuh Hamiltonian of the system. 

Near the point where the bands of graphene touch each 
other (Dirac point) a simphfied Hamiltonian of a low- 
energy quasiparticle moving in a slowly varying static 
potential U (r) and interacting with an external electro- 
magnetic radiation can be written for a single valley and 
for a certain direction of spin aP^ 



n 



va{y--^A{t))^U{v). (3.2) 



Here p is the momentum of the quasiparticle, v- the 
Fermi velocity, and a- the vector of the Pauli matrices 
in the sublattice space ("pseudospin" space). 

The electromagnetic radiation is taken into account 
choosing a proper vector potential A(t). For a linearly 
polarized monochromatic electromagnetic wave it can be 
taken in the form 



A{t) = -Ecos(m), 



(3.3) 



where E is the amplitude of the electric field in the wave. 
Substituting Eq. (3.2) into Eq. (3.3) and using 



Eq. (3.1) we obtain a complete description of electron 
motion of graphene. In the next three subsections we 
study the properties of solutions of the time-dependent 
Schrodinger equation in graphene in presence of linearly 
polarized monochromatic radiation. 



A. Dynamical gap in the quasiparticle spectra 

As it has been shown long ago, in the Ref. 27, reso- 
nant interaction of quasiparticles with EF can lead to the 
formation of a dynamical gap in the electron spectrum of 
semiconductor. The dynamical gap occurs in the vicinity 
of the values of momentum p, determined by the reso- 
nant condition ec(p)— e^(p) = h£t^ where e^(p) and ec(p) 
are the electron energies in the valence and conduction 
bands respectively. 

In order to illustrate this, let us consider an electron 
propagating in the absence of external potential poten- 
tial \U (r) = 0] along a certain z axis with momentum p, 
perpendicular to the field E, which is assumed to be di- 
rected along the x axis in the plane of the graphene sheet 
(Fig.[l}. 

Let us emphasize, that our choice of the coordinate 
system and associated pseudospin basis differs from that 
usually chosen in graphene-related papers. Usually, all 
the calculations are done in the basis of states, defined 
on different graphene sublattices, so that the z axis is 
directed perpendicularly to the graphene plane, while x 
and y belong to the plane. However, for the purposes of 
the present paper it is more convenient to work with the 




FIG. 1: (Color online) Graphene junction in presence of ex- 
ternal radiation. Axis z is directed along the electron momen- 
tum p in the graphene plane, axis x- parallel to the electric 
field E in the external electromagnetic wave. 



Hamiltonian, having a diagonal form in the basis of states 
{"pseudospin parallel to the momentum p", "pseudospin 
antiparallel to p"}. 

For this choice of the coordinate system the Hamilto- 
nian [Eq. (|3.2[)] takes the form 



where 



H = vpz^z + 2Aaa; cos(r2t). 



v\e\E 



2Vt 



(3.4) 



(3.5) 



As it will be clear from the analysis below, even a weak 
EF (A <C hO) strongly affects the motion of electrons 
with momenta close to the resonant values. 



Pr 



2^' 



(3.6) 



which makes this region of momenta the most interest- 
ing. Therefore, we consider now a quasiparticle with mo- 
mentum pz close to one of these resonant values, say, 
to hQ/{2v). Applying the rotating- wave approximation 
(RWA]£^, we decompose the second term of the Hamilto- 
nian ( |3.4[ ) into the right and the left circularly-polarized 
waves, rotating about z axis in the pseudospin space, 
and neglect the second one in the vicinity of the reso- 
nance chosen. In the other words, we make the following 
replacement in Eq. (3.4): 



2Ad-x cos{Qt) — > 
A{d-x cos{Qt) + d-y sm{Qt)). 



(3.7) 



Then, going to the reference frame (RF) rotating coun- 
terclockwise about the z axis with the angular velocity 
Q in the pseudospin space, we come to a problem with 
a static Hamiltonian. This procedure corresponds to the 
application of the unitary transformation 



U 



(3.8) 



to the two-component wavefunctions. In the new basis 
we obtain the time-independent effective Hamiltonian 



n' 



VPz 



(3.9) 
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The eigenvalues of the Hamiltonian (3.9) determine the 
quasiparticle spectrum as 



VPz 



(3.10) 



Equation (3.10) shows that the spectrum acquires a 



dynamical gap 2A at the resonant value of momentum. 
This gap is proportional to the amplitude E of the ra- 
diation and inversely proportional to its frequency Vt. 
Far from the resonance (Ip^'i; — fil7/2| ^ A), the cor- 
responding quasiparticles are just conventional electrons 
and holes in the absence of radiation, having spectra 
-^\vpz — fiVt/2\ in the chosen rotating RF. 

When the momentum pz approaches the other reso- 
nance, —Ml/2v^ one can analogously calculate the quasi- 
particles spectra in the RWA in the RF rotating clockwise 
about the z axis, 



VPz 



my 



(3.11) 



Let us emphasize that the applicability of the RWA 
for the derivation of each of the Eqs. (3.10) and (3.11) 



is limited not by the closeness of the momentum to the 
corresponding resonance, but by the negligibility of the 
influence of the other resonance. For instance, Eq. (3.10) 
is valid as soon as 



Pz 



2^ 



> 



(3.12) 



The dynamical gaps in the quasiparticle spectra have 
been calculated assuming that EE was linearly polarized 
perpendicular to the direction of momentum p. If the 
field E is directed at some angle 7 with respect to its 
transverse component E^, perpendicular to p, then the 
effective dynamical gap decreases and becomes equal to 



A(7) = A cos 7. 



(3.13) 



As we will show in Sec. IVIIl for the reasonable values 



of the radiation power this gap is much smaller than all 
the other energy scales in the problem, such as the Fermi 
energy, the photon energy hQ^ the heights of the potential 
barriers, and the typical voltages applied to a sample. 

Radiation-induced gaps have been observed in the 
spectra of spontaneous radiation of conventional semi- 
conductors subject to a strong monochromatic laser 
field^^ In this paper we consider the effect of radi- 
ation on the transport properties of graphene. Measure- 
ment of the current is a simpler task and we hope that 
the corresponding experiments are also possible. 

In the limit of small value of A considered here (A <C 
hft), the conductivity of a homogeneous irradiated semi- 
conductor sampl^^or the conductance of a tunnel junc- 
tion between two irradiated uniform sample^^ can be 
strongly affected by the radiation only provided the 



Fermi level is close to the position of the EE-induced gap. 
Any considerable spatial variation of potential due to in- 
homogeneities in the system smears the gap and makes 
its observation in bulk semiconductors very difficult. 

In the present paper we consider an essentially different 
situation when electrons move in a non-uniform step-like 
potential and their chemical potentials in the left and the 
right reservoirs are not necessarily close to those, corre- 
sponding to the resonant momenta. Nevertheless, for a 
broad interval of energies of the incident electrons each 
of them achieves resonances in one or several resonant 
points. 

The reason for this unusual behavior is the coordinate 
dependence of the electron momenta that determine the 
behavior of the wavefunctions. The momenta of electrons 
reach their resonant values near the step. The region, 
where these momenta correspond to the motion inside 
the dynamical gap, is thus localized in space close to 
the step-like potential. Electrons must tunnel through 
the gap in order to contribute to the current between 
the leads. In the regions between the resonant points 
electrons weakly interact with the radiation and prop- 
agate freely. However, the current- volt age characteristic 
of such a junction is determined by the tunneling through 
the dynamical gaps. 



B. Dynamics of electrons normally incident on 
irradiated potential barrier 

In this and the next subsections we study the tunneling 
through the dynamical gap in the momentum space in the 
vicinity of the resonant point. 

Consider first an electron normally incident on a 
smooth potential barrier U{r) = U{z), varying only in 
one dimension. To be specific, we assume that the elec- 
tron is moving in the conduction band along the z axis, 
and the height of the potential barrier increases with z. 
Let e be the total electron energy (kinetic+potential) in 
the initial laboratory RF far from the resonant point zq 
[determined by the equation U{zo) = e — where the 
momentum in the absence of the radiation would equal 
to the resonant value hQ/2v (see Fig. [2|. The EE in the 
graphene plane is polarized perpendicularly to the z axis, 
the direction of electron momentum. 

As follows from the Hamiltonian ( |3.2[ ), the electron 
velocity far from the resonance, where \U{z) — U{zo)\ 
A, is determined by its pseudospin: 



V = va. 



(3.14) 



This means that far from the resonant point the veloc- 
ities and the pseudospin | |) of the incident e*^ f^p{z)dz 
and the transmitted te*^ p{z)dz ^g^y^g ^-^^ both di- 
rected along the z-axis. The momentum p{z) entering 
the exponents is determined by the equation 



p{z) = {s-U{z))/v. 



(3.15) 



6 




^ij p(z)dz ^ 

incoming wave 



reflected 



zo 

resonant point 



FIG. 2: Scattering of quasiparticles moving perpendicularly 
to the barrier. 



Let us emphasize, that in contrast to the correspond- 
ing quasiclassical expression for the Schrodinger equa- 
tion, Eq. (3.15) is exact, i.e. valid for arbitrary values of 
momenta p, which follows immediately from the Dirac- 
type Hamiltonian (|3.2[). According to our notations, the 



pseudospin | i) of the reflected wave is antiparallel to the 
z axis. 

Let us find the transmission and reflection coefficients 
of the electron in the region of the barrier where it 
strongly interacts with the radiation. In the pseudospin 
space we go to the rotating RF defined by transformation 
(3.8). In this frame the Hamiltonian is static, the particle 



has energy e — ^ and scatters elastically. The transmis- 
sion and reflection coefficients are found solving the Dirac 
equation for the quasiparticle wavefunctions close to the 
resonant point. In the presence of the potential U {z) this 
equation takes the form 



VPz 



^ = 6^. (3.16) 



Here ^ is the two-component wavefunction. 

Without losing generality, we can set zq = 0. Since the 
potential is smooth, it can be expanded in small z and 
become linear around the resonant point: 



U{z) 



-Fz. 



(3.17) 



The Dirac equation for the particle with energy e - 
can be written in the momentum representation as 



nn 

2 



op 



Aa, ^(p)(3.18) 



We have taken into account that in the momentum rep- 
resentation z = ih-§-. 

op 

Equation (3.18) describes Landau-Zener tunneling 
through the dynamical gap (Fig. [3| in the momentum 
space. This phenomenon is analogous to the electron 
tunneling through the forbidden band in a conventional 
semiconductor tunnel p-n junctions^^ or through a non- 
irradiated graphene p-n junctions, when the incident elec- 
tron has some finite transverse component of momentum. 




FIG. 3: Tunneling through the dynamical gap in the mo- 
mentum space (rotating R F). T he quasiparticle spectrum in 
presence of radiation, Eq. (3.10), is shown by the solid lines 
in the plot, the spectrum in absence of radiation- by the thin 
lines. 



The latter case, when the role of the gap is played by the 
quasiparticle energy at zero longitudinal momentum, has 
been studied in the Refs. 41 and 42, 



Equation (3.18) can be solved exactly and the tunnel- 
ing probability T = |tp takes the form 



■ hvF 



(3.19) 



We emphasize that Eq. (3.19) gives the exact probabil 



ity of Landau-Zener tunneling, i.e. it is valid for arbitrary 
values of the parameter 



C = 



hvF 



(3.20) 



in the exponent rather than only in the "quasiclassical" 
limit (£ ^ 1), considered in the previous paper J^^^. 

According to the energy conservation law in the rotat- 
ing RF, the reflected wave has the momentum ^ — p{z) 
at the same point z^ where the incident wave had mo- 
mentum p{z) (Fig. |3|. This momentum corresponds to 
the energy e — HQ of the scattered electron in the ini- 
tial laboratory RF. This means, that reflecting from the 
resonant point the electron emits a photon of the energy 

nn. 

Similarly, one can consider the cases of the quasipar- 
ticle scattering in the valence band of graphene at the 
resonant points, where the momentum in absence of ra- 
diation equals the other resonant value — 

So, we conclude that an electron normally incident on 
an irradiated potential barrier can pass with the prob- 
ability T through the resonant point, where its kinetic 
energy vp in absence of radiation would be equal HQ/ 2 or 
be reflected from that point with the probability 1 — T, 
emitting a photon of energy hft. At the other type of 
the resonant points, where the kinetic energy in the ab- 
sence of radiation equals the particle absorbs a 
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photon of energy HQ during the reflection. The emis- 
sion/absorption is accompanied by the particle reflection 
and the pseudospin-flip. 

Tunneling suppression. Applying strong radiation or 
smooth potential barriers, such that C 1, o ne ca n 
achieve small tunneling probabilities T <C 1 [Eq. ( |3.19 )]. 
Such a situation was considered in the Ref. 36, where it 
was argued that in this case one would be able to sup- 
press Klein tunneling and conflne electrons by irradiated 
barriers. In Sec. IVIII we revisit this issue and estimate 



the radiation power needed to realize such a conflnement. 

Limit of small C. Let us show now, that in presence of 
impurities, elect ron-phonon and electron-electron inter- 
actions in the sample in the limit £ <C 1 the probabilities 
of reflection and transmission at the resonant point re- 
main the same as in the ballistic sample even for the 
radiation intensities S much below the critical value Sc^ 
at which the gap is suppressed in a wide uniform sam- 
ple. Making this statement we imply that C is given 
in Eq. (3.20) with A being not the real dynamical gap. 



but the quantity deflned by Eq. (3.5), which gives the 



value of the gap for a clean sample neglecting impurities, 
electron-phonon and electron-electron interactions. 

For small jC Landau-Zener tunneling in the ballistic 
sample occurs in the momentum interval of the order of 



A Plz 



(3.21) 



In the coordinate space the corresponding tunneling 
length is 



A vlz 



V A Plz 



(3.22) 



Note, that the scales (3.21) and ( 3.22[ ) are determined by 
the potential slope only and do not depend on the value 
of the gap. 



Equation (3.22) shows that the perturbation A&x in 



the quasiparticles Hamiltonian ( |3.4| ) induces transitions 
between the states | j) and | |) in a wide spatial interval 
independent of A. This occurs even when the dynami- 
cal gap in the particle spectrum is smeared due to the 
radiation-independent relaxation processes. 

The effect of the radiation-independent relaxation on 
the Landau-Zener tunneling can be neglected provided 
the probability of excitation due to it is much smaller 
than 1, or, in other words, the relaxation time tr is suf- 
flciently large 



TR > 



(3.23) 



This condition does not necessarily imply that the prob- 
ability of reflection C = ii/S? / {hvF) in the Landau- 
Zener tunneling is much larger than the probability of 
an EF-independent transition on the corresponding in- 
terval. However, using the inequality (3.23) will allow us 



to consider the EF-independent transitions and the tun- 
neling independently of each other in the calculation of 
the photocurrent in the subsequent sections. 



C. Radiation-induced hops between trajectories at 
the resonant points 

In the previous subsection we considered electrons in- 
cident normally on a potential barrier in a transverse EE. 
In this subsection we analyze the scattering in the pres- 
ence of radiation for arbitrary angles between the elec- 
tron momentum p at the resonant point, the slope of the 
potential dU/dr^ and the fleld E. 



As it has been already mentioned [cf. Eq. (3.13)], the 
effective dynamical gap close to the resonant point is A = 
A cos 7, where 
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2 



(P,E) 



(3.24) 



is the angle between the field E and its transverse com- 
ponent E^, perpendicular to the momentum p. 

The effective potential slope along the direction of 
the electron momentum at the resonant point is F = 
I F cos/? I, where 



[3 = {dU/dr, p) 



(3.25) 



is the angle between the slope of the potential and the 
momentum at the resonant point. 

Then, instead of Eq. (3.19), one should use the formula 



g hv\Fcos(3\ 



(3.26) 



for the probability for an electron to pass through the 
resonant area without being scattered by the radiation. 
Below we will present an explicit derivation of this result. 

Note also that, although being reflected at the reso- 
nant point the electron changes its pseudospin (direction 
of velocity), it does not necessarily follow the same tra- 
jectory along which it moved towards the resonant point. 
The reason is that the reflection is accompanied by a pho- 
ton emission or absorption, and the new electron energy 
^refi = s ^ HQ corresponds to a new trajectory. This 
means that the radiation can enforce the electron to hop 
onto the trajectory, corresponding to the opposite pseu- 
dospin at the resonant point (Fig.[4|. Let us now proceed 
to an explicit derivation of this result and of Eq. (3.26). 



As discussed in the previous subsection, far from the 
resonances the electron motion can be considered qua- 
siclassically neglecting the radiation. Let AOB (Fig. [4| 
be a classical trajectory of a quasiparticle moving in ab- 
sence of radiation in the potential U{r). Here, O is the 
resonant point where the momentum equals the resonant 
value p = hQ/{2v). At the resonant point the particle can 
be reflected and hop to the other trajectory COD. Tra- 
jectory COD corresponds to the same momentum but 
the opposite pseudospi n at the resonant point (Fig. |4|. 

To obtain the Eq. (3.26) and analyze the radiation- 



induced hops between trajectories we introduce at each 
point of AOB the coordinate system r^C^, where the axis 
^ is directed tangentially along the trajectory, (- nor- 
mally to the graphene sheet, and 77 lies in the plane of 
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FIG. 4: (Color online) Radiation-induced hop to the trajec- 
tory with opposite pseudospin at the resonant point. Moving 
initially along the arc AO, at the resonant point O the quasi- 
particle either hops to the trajectory OD, if scattered by the 
radiation, or continues moving along OB, in case the scat- 
tering has not occurred. Axis ^ of the coordinate system is 
directed along the trajectory AOB at each point, axis ^ is 
perpendicular to the graphene plane, and 77 lies in the plane. 



the sheet. Let s be the particle coordinate along the tra- 
jectory AOB. Introduce also the fixed coordinate system 
xyz, such that the x axis is directed along the EF and y- 
normally to the graphene plane and along (. The system 
r]((, is obtained rotating xyz-system around the y axis by 
some angle a{s). 

The corresponding pseudospin transformation is 



(3.27) 



Let us set 5 = at the resonant point. The Hamito- 



nian (3.2), written in the pseudospin basis defined by the 



system r]((_ and linearized in the vicinity of the point O, 
takes the form 



dU 



dU 

+2A(a2 sin a + cFx cos a) (zo^(Vtt). 



(3.28) 



The third line of Eq. (3.28) describes the motion of the 



wavepacket in the direction of axis 77, perpendicular to 
the trajectory. This part of the Hamiltonian is indepen- 
dent of the longitudinal coordinates on the scales being 
considered and can be disregarded if we want to find the 
probability of tunneling through the point O along the 

trajectory AOB. 

The term —h{da/ds)ay in the second line of Eq. (3.28) 



is the gauge potential generated by the local rotations of 
the reference frame 77^^. Considering the tunneling in 
the vicinity of the point O, one can neglect this term, 
provided that the potential is smooth enough {X/v 
r^, where A is the characteristic scale of the potential 
variation). Note, that such a neglect in general can be 
made only on a small enough interval, but not on the 
whole trajectory between different resonant points. 

In the related paper, Ref. |36l the local reference frame 
reverted the direction of two axes at the point where p = 



0, but the corresponding non- negligible delta-function- 
type gauge potential, generated by this rotation, was not 
taken into account. 



The last line of Eq. (3.28) is the EF-induced time- 



dependent perturbation of the Hamiltonian. Following 
the same line of reasoning as when deriving RWA^^, one 
can keep in this perturbation only the transverse circu- 
larly polarized wave rotating clockwise about p (axis ^). 

Making all the aforementioned approximations in the 
Eq. (3.28) and going to the rotating RF, we arrive at the 
same problem as considered in the previous subsection 
with the effective gap A cos 7 and the eff ective potential 
slope Fcos/3, and finally obtain the Eq. (3.26). 



IV. DYNAMICS OF ELECTRONS IN A 
BALLISTIC GRAPHENE P-N JUNCTION 

From now on we restrict ourselves to the consideration 
of a specific type of potential barriers. We assume that 
the potential U{z), varying only along a certain axis 
increases monotonically from U{—oo) = far in the left 
lead to U{-\-oo) = Uq > far in the right one. The 
electric field E of the external radiation is directed along 
the X axis, perpendicular to the z axis and lying in the 
graphene plain. 

From the previous section we know that the electron 
motion in a non-uniform potential in the presence of ra- 
diation can be considered as quasiclassical between the 
resonant points, where the hops between trajectories can 
take place, accompanied by the pseudospin-flips and pho- 
ton emissions or absorptions. 

Now we have to consider all the classical electron tra- 
jectories in the potential under consideration and find 
their resonant points. 



A. Dynamics of electrons non-interacting with 
radiation 



In the absence of radiation electron transmission 
through the junction is determined by its transverse mo- 
mentum pj_ = psin^, which is conserved during the mo- 
tion. The electron tunnels through the non-irradiated 
p-n junction with the probabilit}^ 



To = e- 



(4.1) 



and, accordingly, reflects back with the probability 1 — Tq. 
Here Fq is the effective potential slope at the p-n inter- 
face, calculated taking into account the charge density 
distribution in the junctioiP^. For each incident electron 
Fo should be understood as the slope of the potential at 
the point where the longitudinal momentum equals zero. 
Note, that in general it is different from the potential 
slope F at the resonant point. 

Let us classify the regimes of the particle motion 
in the absence of radiation according to their trans- 
verse momenta. We introduce the ^'2D-modes'\ with 



9 



\p±\ ^ {^Fq / {ttv)) ^ ^ and the ''normal modes^\ with 
\pi_\ <C (fiFo/(7r'u)) 2 . In the first case, for large enough 
p^, the electrons perfectly reflect from the interface, 
from the place where their longitudinal momentum turns 
to zero. In the second case particles freely penetrate 
through the junction without reflection. 

The classical path of an electron in the presence of 
the radiation consists of the pieces of its paths in the ab- 
sence of the radiation, stuck to each other at the resonant 
points. 

Let us show now that on each electron trajectory in the 
presence of the radiation there can be no more than two 
resonant points. An electron encounters a resonant point, 
when its full momentum (p^ + ) ^ reaches the resonant 
value Ml/(2v)^ i.e. when the longitudinal momentum 
becomes equal to one of the two values 



= ± 



(4.2) 



The momentum pz of a classical electron, moving in the 
monotonically increasing potential U{z)^ changes in time 
monotonically, Pz = —dzU < 0. If the particle scatters at 
a resonant point, its momentum does not change. Then 
each of the momentum values (4.2) of Pz is reached not 



more than once, and the electron encounters correspond- 
ingly not more than two resonant points. 

In order to calculate the total current through the junc- 
tion we have to consider the electron paths starting from 
the left lead and terminating in the right lead. The con- 
tribution to the current of the inverted processes, trans- 
mission from right to left, has been already taken into 
account during the derivation of Eq. (2.7). 



B. 



Electron paths in effectively two-dimensional 
modes. 



Consider, for insta nce, the process illustrated in the 
Fig. [sja). Spectrum (4.3) is shown by the dashed lines 
there. The electron path on this plot, indicated by solid 
lines with arrows, starts at the branch of the spectrum 
corresponding to the conduction band and positive veloc- 
ity Vz = de/dpz' If there was no radiation, the electron 
would always be in the conduction band and would follow 
the dotted line in the figure. In the end of the dotted line 
Vz < 0^ indicating the fact that the electron would have 
reflected from the junction. In fact, when reaching the 
first resonant point, in the process under consideration 
the electron emits a photon [curved line in the Fig.jsja)], 
and then moves in the valence band, passing through 
the second resonant point. In the end of the process the 
sign of the longitudinal velocity coincides with the ini- 
tial one. This means, that as a result of the process, 
the electron penetrated through the junction from one 
lead into another. Thus, the radiation can assist electron 
transmission through the junction. 




Pres 




b) 



FIG. 5: (Color online) Process contributing to the current 
through the junction (solid lines). The particle emits a pho- 
ton at the first resonant point, and passes without reflection 
the second ("2D-modes"). (a) In the plot "kinetic energy vs. 
longitudinal momentum" . (b) In the spatial coordinate space 
(graphene plane view) . The emission of the photon is accom- 
panied by the hop at the resonant point between two classical 
paths, obtained in absence of radiation. 



Let us find first the resonant points for an electron 
in the "2D-mode", i.e. having large \p±\, and incident 
on the barrier from the left with a positive energy and 
momentum p > hft/{2v). 

If \p±\ > hft/{2v), then there are no resonant points 
on the trajectory. The electron weakly interacts with EF, 
perfectly rebounds from the p-n interface in presence of 
radiation and thus does not contribute to the current 
through the junction. 

If \p±\ < hVt/{2v) and the potential is high enough, 
then there are two resonant points, corresponding to the 
values (4.2) of momentum p^. 



The behavior of the particle is best illustrated as its 
path on the plot of its kinetic energy in the laboratory 
RF, 



^kin ijPz 



2^i 



-Pi) 



(4.3) 



versus the longitudinal momentum Pz . 



Let us now take a look at the corresponding electron 
trajectory in the coordinate space. Fig. jsjb). Shown 
there two classical trajectories in absence of radiation 
are tangent at the resonant point. The left trajectory in 
the absence of radiation starts and ends up in the region 
of n-type graphene. It is the path of the electron in the 
conduction band, having classical Hamiltonian 



^con(i(p,r) =vp^U{y). 



(4.4) 



Similarly, the right path begins and terminates in the p- 
type graphene, and corresponds to the motion of a par- 
ticle in the valence band with the Hamiltonian 



7Y^az(p,r) = -vp^U{y). 



(4.5) 



Due to the effect of the radiation the particle hops at the 
resonant point from the left path to the right one, the 
whole resulting trajectory is shown by the solid line in 
the Fig.jsjb). 
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For a given electron energy and the angle of incidence 
the probability to pass along the trajectory under con- 
sideration is 



Plscatt,2pass — (1 ~ ^l)^2 5 



(4.6) 



where Ti and T2 are the transmission probabilities cor- 
respondingly at the first and the second resonant points. 

Analogous to Fig. [Sj Fig. [6] illustrates the tunneling, 
when the particle is not scattered at the first resonant 
point, but rebounds from the second one, emitting a pho- 
ton. 




'^1 




resonant 







s^_^^---7 points 


-type 




■typei 









b) 



FIG. 6: (Color online) Another process contributing to the 
current through the junction (solid lines). The particle passes 
without reflection the first resonant point and emits a photon 
at the second ("2D-modes"). (a) In the plot "kinetic energy 
vs. longitudinal momentum", (b) In the spatial coordinate 
space (graphene plane view). The emission of the photon is 
accompanied by the hop at the resonant point between two 
classical paths, obtained in absence of radiation. 

Let us emphasize, that one should not confound the 
static gap 2v|px| of the Hamiltonian 



(4.7) 



giving the particle spectrum in the absence of radiation 
in the Figs. [5] and [6j with the dynamical gap 2A in 
Fig. |3) where the quasiparticles spectrum is shown in 
the rotating reference frame, taking into account the ra- 
diation. Landau-Zener tunneling through the static gap 
is strongly suppressed, since \p±_ \ ^ {HFq/ (ttv))^ . 

In the processes shown in the plots in Figs, [sja) and 
[6]^a), electron ends up in the right lead, because finally 
the particle is on the branch of the spectrum correspond- 
ing to the same direction of velocity at \pz\ -^00 as on 
the initial one. In contrast to that. Fig. [7| illustrates the 
backscattering: the particles enter and leave the junction 
in the same lead L. According to the results of Sec. [ll| 
the corresponding trajectories do not contribute to the 
current and should be excluded from our consideration. 

We see that, although in the absence of radiation elec- 
trons in the normal modes could not tunnel through the 
junction, they can penetrate from the left to the right 
lead emitting one photon when the sample is irradiated 
by EF. Due to the time-reversal symmetry, the inverse 
processes, when a particle penetrates from the right to 
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a) 



b) 



FIG. 7: Kinetic energies of the particles incident from the 
left reservoir and scattered back there ("2D-modes"). a) The 
particle emits and absorbs a photon at the first and the sec- 
ond resonant points correspondingly. b)The particle is not 
scattered by the radiation. 



the left lead absorbing a photon, also exist. As we dis- 
cussed before, in order to derive the current- voltage char- 
acteristics of the junction, we have to calculate the tun- 
neling probabilities only for the particles incident from 
the left, the contribution of inversed processes has been 
taken into account when deriving Eq. (2.7). 



C. Electron paths in normal modes. 

Now let us consider the normal modes, i.e. having such 
small transverse momenta that in the momentum space 
the particles freely penetrate through the static gap, as 
if it was zero. This corresponds to p± = in Eq. ( |4.3[ ). 

For the tunneling from the left to the right lead both 
possible electron paths €kin{Pz) are shown in Fig. [s] 



Ski', 




Ski', 
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a) 



b) 



FIG. 8: Kinetic energies of the particles incident from the 
left reservoir and penetrating into the right one (normal 
modes). a)Usual Klein tunneling, no radiation-induced scat- 
tering, b) Tunneling accompanied by the two-photon emis- 
sion. 

The first one, in the Fig.jsja), corresponds to the usual 
Klein tunneling, when the particle is not scattered by the 
radiation. The static gap is negligible, so the particle 
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passes from the conduction to the valence band through 
the Dirac point without reflection. 

The second path, in the Fig. |8|b), illustrates the pro- 
cess of tunneling from the left to the right lead, accompa- 
nied by two photon emissions. Note, that since the parti- 
cle is being twice reflected by the EF in this process, the 
radiation, whatever strong, cannot suppress this channel 
of tunneling. The increase in the radiation power leads 
only to the enhancement of the transmission probability 
in the channel. 



V. PHOTOCURRENT IN BALLISTIC SAMPLES 

A. Photocurrent due to effectively two-dimensional 
modes 

For the experimentally relevant parameters (see Sec. 



according to Eq. (2.7), only the electrons with energies e 



VII for more details), the electron Fermi momentum in 



the leads is much larger than the characteristic transverse 
momentum needed for a particle to be reflected at the p-n 
interface, 



Pf > 



TTV 



(5.1) 



Hence, the majority of the particles are in the "2D- 
modes". In this subsection we calculate the photocur- 
rent in a ballistic graphene p-n junction, neglecting the 
normal modes. 

As we have just shown for such modes, the tunneling 
from left to right in a monotonically increasing potential 
is necessarily accompanied by the emission of one photon. 
Since the tunneling from left to right involves photon 
emission, tunneling from right to left, due to the time- 
reversal symmetry, involves photon absorption. 

Then for some energies 

PLR{e, e - hn, 9) ^ PLR{e - hQ, e, 0) = 0, (5.2) 

and, according to the results of Section [TTj the photocur- 
rent flows through the junction in absence of any voltage 
applied to it. 

If the height of the potential obeys the inequality Uq > 
then at least one of the processes in the Fig. [s] or 
Fig. |6] or time-reversed processes is possible in the p-n 
junction, and some photocurrent flows through it. From 
now on we assume, however, that the barrier height and 
the Fermi energy are large enough. 



Uq - Ef > 



(5.3) 
(5.4) 



As it will be clear from the further consideration, under 
these conditions the photocurrent is maximal. 

Since we are neglecting normal modes, and in the "2D- 
modes" the tunneling from left to right is possible only 
with a single-photon absorption, in Eq. (2.7) we have to 
discard all the terms except those with n = — 1. Then, 



in the left lead, such that £ — hft<eF<£,OT^ the same, 
Sf < £ < Sf contribute to the current. 

This can be understood as follows. Each electron with 
energy sf — hfl < e' < Ef '^^ the right lead penetrates 
into the left lead increasing its energy by HQ (Fig. [9|. 
The current, carried by these electrons, is not compen- 
sated by the corresponding time-reversed processes, since 
the states of electrons, incoming from the left lead on 
the energies Ef < £ < £f + ^^ above the Fermi level, 
are unoccupied. So, the electrons coming from the leads 
on the correspondingly lower energies either compensate 
each others' contributions to the currents or return back 
to their initial leads. 



n 




R 



FIG. 9: (Color online) Contribution of electrons with dif- 
ferent energies into the photocurrent through graphene p-n 
junction. 

Thus, in order to fin d the photocurrent in the junction, 
one has to apply Eq. (2.7), performing integration over 
the states of the left lead with energies in the interval 
{sf;£f + ^^). 

Let us assume that at some place around the resonant 
points the potential has a constant slope F. Then, ac- 
cording to the results of Section |III[ the probability for 
an electron at some energy £ and angle of incidence to 
go along the trajectory under consideration, is 

= T(l - T) = e-^^°^^ (1 - e-^^°^^) , (5.5) 

where 7 is the angle between the electron momentum at 
the resonant point and the normal to the junction [cf. 
Figs. [5|b) and[6]^b)], C is the parameter introduced by 
Eq. (IsM. 



According to Eq. (2.7) the photocurrent is 



^2D 



{£F-\-hn)/v 



dp 



£f/v 



2pv 



pdO 
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The angle 7 can be determined from the transverse mo- 
mentum conservation law 



psmt 



2^ 



sm7. 



(5.7) 



Then 



I2D = ^eW / dp (i7 ^^_^^^ cos 7 

Jef/v J — 7t/2 



xe 



-C cos 7 



(1 



-C cos 7\ 



(5.E 



■KV 



{Li{C) - h{C) - ii(2£) + Ji(2£)) , (5.9) 



where Li{z) and are correspondingly the modi- 

fied Struve and the modified Bessel functions of the first 
order^^. 



According to the Eqs. ( |3.2Q[ ) and ( |3.5[ ), the parameter 

(5.10) 



£ is proportional to the intensity S of the EF: 
^^e'v S 



hcF n'^' 



Expression (5.9) is plotted in Fig. 10 
in = . 



(5.11) 



The probability to tunnel in each channel, i.e., along each 
trajectory in the Figs.jsjb) and[6]^b), is £, in the limit of 

small C. Then the Landauer-type conductance of each 

2 

channel is G = Since each electron acquires the 

energy HQ tunneling from right to left, the effect of ra- 
diation on the conducting channels is equivalent to the 
effective voltage Veff = hQ/e applied to the junction. 
Then Eq. (5.12) transforms into 



TsmaLi 
^2D 



N2DGVeff. 



(5.14) 



Large-intensity regime. 



In the limit of large radiation powers (£ ^ 1) Eq. (5.9) 
reduces to 



jlarge 



ZeWQ? 1 



(5.15) 



The photocurrent due to the modes under consideration 
is being suppressed by strong radiation, since the elec- 
trons tend to be refiected at each resonant point when 
C is large and the tunneling along the paths in the 
Figs.jsjb) and[6]^b) becomes unlikely. 

Since the photocurrent carried by electrons with suffi- 
ciently large transverse momenta vanishes at £ ^ 00, one 
needs to consider the contribution of the other electrons 
incident almost normally at the p-n interface. 



0.05 -- 




FIG. 10: Photocurrent in a graphene p-n junction as a func- 
tion of the radiation intensity. 

Small-intensity regime. 

In the experimentally relevant regime of small radia- 
tion intensities (£ <C 1) the photocurrent is 



Tsmall 
'-2D 



27TV 



C. 



(5.12) 



The last formula can be understood as follows. The effec- 
tive number of conducting "2D- modes" , i.e. of transverse 
channels in the energy interval ~ H^l, is 



N2D = W 



(5.13) 



B. Photocurrent due to normal modes 

In the normal modes electrons have very small trans- 
verse momenta p_\_, so the problem can be viewed as one- 
dimensional. Introduce the effective number of normal 
modes, i.e. of the channels where electrons propagate 
without reflection at the p-n interface (taking into ac- 
count spin and valley degeneracies): 



27rh 



dp± e 



^0 



2W lhF( 



V 



. (5.16) 



As discussed before, when the condition (5.1) is satisfled, 
one has the inequality N2D ^ ^iD- 

The corresponding effective one-dimensional density of 
states (per a unit of longitudinal length) is 



^"1 



ID 



(5.17) 



The current, flowing through the juncti on w ithin the 
one-dimensional picture, i^^ [cf. also Eq. (2.7)] 



/= / de viD PLi^(g, g + nhVt) 

^ n 



(5.18) 



where Plr{£^£ + nhQ) is the probability to tunnel from 
left to right from the state with energy £ to the one with 
e + nhft. 
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The only tunneling process contributing to the pho- 
tocurrent [Fig. [sj^b)] is the one when the particle is being 
reflected at both resonant points, emitting two photons. 
The probability of this double reflection is (1 — e~^)^. 

The photocurrent due to the normal modes is possible 
provided Uq > hVt. Again, we assume that the height of 
the potential barrier and the Fermi level are large enough, 
so that the photocurrent due to the processes under con- 
sideration not simply exists, but also reaches maximum 
as a function of these parameters. 



Sf > 



Uq-£f> 



3hn 
3hn 



2 ' 



(5.19) 
(5.20) 



so that both resonant points are present on the potential 
for each electron in the energy interval of the width 2hQ. 
Then the photocurrent due to the normal modes is 



_£.2 f2hQ\ 
e 



(5.21) 



We see thus, that at very large intensities {C 1) the 
photocurrent is saturated. This happens because in this 
limit any electron, reflected twice by the radiation with 
the probability close to 1, recovers its initial direction of 
velocity and penetrates into the opposite lead. 



C. Full current vs. radiation intensity 

To sum up, the full photocurrent / flowing through the 
irradiated p-n junction is given by the sum 



I = L 



2D 



ho 



(5.22) 



with I2D and Iio given by Eqs. (5.9) and (5.21) respec- 
tively. 



If the radiation power is not too large. 



3 



2D 



47r NiD ' 



(5.23) 



the p-n junction is effectively two-dimensional, and the 
photocurrent in it is described by Eq. (5.9). For larger in- 



tensities, when the last inequality reverts, the photocur- 
rent strongly decreases (by the factor of ~ N2d/Nid) 
and saturates at the value Iid{^ 00) = ehQNuj/TT. 



D. Suppression of the tunneling 

Current-voltage characteristic of a strongly irradiated 
junction. Assume now that the voltage V is applied to 
the junction. Here we make the convention that the volt- 
age is positive in the case of the forward bias, when the 
electric potential of the p-type graphene is larger than 
that of the n-type. The corresponding direction of the 
current in the absence of radiation is considered positive. 



Introducing the conductance Gi^aii of the ballistic 
graphene p-n junction 



Gba 



h 



(5.24) 



and taking into account Eq. (5.21), the current- volt age 



characteristic of the junction under a very strong irradi- 
ation (£^ N2d/Nid) can be written as 



liV) = Call V 



(5.25) 



The second term in the parenthesis of the last equation is 
the photocurrent, Eq. (5.21) written in the limit £ ^ 00. 
The equation holds until V < Uq— Sp For larger 

voltages the current grows even more rapidly due to the 
appearance of normal modes without resonant points. 

Thus, the current through the p-n junction with high 
enough potential barrier cannot be suppressed by the ra- 
diation, due to the existence of the channels where elec- 
trons can tunnel with probability 1, gaining or losing two 
quanta hVt. As a result, the differential conductance of 
the junction is the same as in the absence of radiation. 

Conditions for tunneling suppression. 

In order to suppress the tunneling in the junction one 
has to cancel the photocurrent due to the "ID modes", 
i.e. to make a sufficiently low potential barrier Uq < HQ 
(or increase the frequency), such that the normal trajec- 
tories with two resonant points no longer exist. 

Making a p-n junction with Uq > hQ/2 and with one 
resonant point for electrons coming with the Fermi en- 
ergy, and then applying strong radiation one can ex- 
ponentially suppress the conductance due to the nor- 
mal modes and suppress the photocurrent carried by the 
"2D modes", because it becomes proportional to oc >C~^, 
Eq. (5.15). One may also think of the situation for 
Uq < hQ/2 when one resonant point still exists for elec- 
trons coming at the Fermi energy. In that case, applying 
a small voltage V to the junction, one obtains the normal 
current / - l[^^^\eV/hQ) - {e^WQ.V /vh)C-^ , which is 
also suppressed at large intensities as ex . 



VI. PHOTOCURRENT IN DISORDERED 
SAMPLES 

In the previous sections we considered purely ballistic 
quasiparticle transport in irradiated graphene samples. 
In principle, the transport can be strongly affected by 
disorder, electron-phonon and electron-electron interac- 
tions. According to the Ref. |46l as concerns their con- 
ductance, junctions fabricated in recent experiments are 
at the best in the crossover between the ballistic and dif- 
fusive regimes. 

In this section we calculate the photocurrent in a p-n 
junction under a weak irradiation, £ <C 1, in a disordered 
sample in presence of electron-electron interaction. We 
show that the photocurrent is the same as previously 
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calculated in a ballistic sample, Eq. (5.12), provided that 
the impurities-induced resistance of some resonant region 
in the junction, where the carriers are effectively excited 
by the radiation, is not large compared to the ballistic 
resistance of the junction. 

In this section we use the following assumptions. We 
neglect the relaxation due to electron-phonon interaction 
on the length of the junction and assume that the con- 
dition (3.23) is fulfilled. The validity of such approxima- 



tions will be confirmed in Sec. |VII| by explicit estimates. 
For simplicity, we suppose that the intervalley scatter- 
ing is weak and neglect it. Taking this scattering into 
account is possible in a similar way, which leads to the 
similar formulas. When writing the kinetic equations we 
also neglect the possibility of disorder-assisted photon 
emission or absorption. In fact, electrons can rebound 
from the impurities, gaining or losing energy quanta HQ 
in the same way as from the smooth potential U (r), con- 
sidered in Sec. Ill Such a process would open additional 
channels for the tunneling of electrons from one lead into 
the another and thus would increase the photocurrent. 

Introduce the distribution function of electrons in the 
conduction band, /|(p,r, t), and in the valence band, 
/l(p,r, t). The first ones have their momenta paral- 
lel to their pseudospins, while the second- antiparallel. 
The stationary distributions are described by the kinetic 
equations (for the derivation see the Appendix) 



.d^ud^ + y^drj /t = r(p)(/^ - /t) + (st/)T, (6.1) 
(^-d^ud^ - v^d„^ fy = r(p)(/t - /i) + (st/)i,(6.2) 

where r(p) is the rate of the radiation-induced 
pseudospin-fiips, and 



St/ = (St/) 



imp 



is the collision integral for 
electron-electron interaction. 



V {Stfr (6.3) 
impurity scattering and 



Kinetic equations (6.1), (6.2) are valid in the case of 
weak enough radiation, £ <C 1. The effect of the radi- 
ation and the effect of other processes on the electron 
dynamics are described by two independent terms in the 
right-hand sides (rhhs) of the equations. The mutual in- 
fiuence of the two corresponding transition rates on each 
other c an be neglected due to the fulfillment of the con- 
dition (3.23), derived in Sec. 



Ill 



from the considerations 
of Landau- Zener tunneling in the momentum space. The 
same condition, 



A24 



follows from the explicit deriva- 
tion of the kinetic equations (cf. Appendix). Let us em- 
phasize, that treating perturbatively the effect of disorder 
and radiation-induced transitions in Eqs. (6.1) and (6.2) 
one should use the values of velocity v and collisional 
terms, renormalized by disordei^^ and electron-electron 
interact ion 

The contribution of elastic impurities into the collision 
integral has the form 

(St/u(p,r))^"'^ 



/ 



rfpVp'(/T,i(p',r) -/T,i(P,r))5(p-p'), (6-4) 



where the quantity Wpp' satisfies the condition following 
from the time-reversal symmetry 



(6.5) 



similarly to the analogous Eq. ([2^ in Sec. [llj 

Analogously one can write the collision integral (St/)^^ 
for the electron-electron interaction (see, for instance, 
Ref. |48|) . However, further we will not need an explicit 
form of this integral using instead only the fact that the 
electron-electron interaction does not considerably mod- 
ify the resistance of the junction, i.e. that such a mod- 
ification is much smaller than the ballistic resistance of 
the junction. 

W e deriv e the photocurrent from the kinetic equations 
(6.1), ( |6.2| ), making perturbation theory in r(p). In the 
zeroth order, when r(p) = 0, neglecting the change of 
the distribution functions due to the collision integral on 
the energy scal es of order of ef and hQ^ we can write the 
solution of Eq. (6.1 ) as an equilibrium Fermi distribution 
/^(p,r,t) = 0{eF — vp — /7(r)), where is the theta- 



function. Similarly, the solution of Eq. (6.2) in absence 
of radiation is (p, r, t) = 6{eF -\- vp — U (r)). 

Taking into account the processes of the first order in 
r(p), one arrives at a slightly modified version of these 
distributions; the radiation excites some small number of 
electrons with energies £F — hVt<£'<£F into the energy 
interval Si? < e < SF^fiQ.. This happens sufficiently close 
to the p-n interface, at \U{z)—eF\ < h£t/2 (see Fig. 11). 
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FIG. 11: Radiation-induced excitation processes, which con- 
tribute to the photocurrent in the graphene p-n junction. 

Since above the Fermi energy the distribution function 
/ is small, / <C 1, and below- / 1, we can rewrite 
Eq. (6.1) for the electrons in the valence band in the 



aforementioned region as 



p 



/^=r(p) + (St/)^ (6.6) 



Let us multiply this equation by the factor 
ApededO / (h'^v) and integrate over the angle ^ or, in other 
words, over the direction of momentum p. Introduce 



s{hn-2pv), 



(6.7) 
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which is the rate r(p), averaged over the direction of p, 
and vq = ^}/{7rhv'^), the density of states per unit square 
at the resonant point (the latter takes into account spin 
and vahey degeneracies). From Eq. (6.6) we get after the 
integration 



de 



eTuQ de -\- e 



d dn^ 



where 



and 



dn^{e,z) 
de 



de 



de = de • 4 



de = de ■ 4e 



pdO 



f^{t,p,z) (6.9) 



pdO 



1^ 

J (2^ 



(6.10) 



are correspondingly the density of electrons and the den- 
sity of the current, carried by electrons with energies in 
the interval {e,e -\- de) in the conducti on b and at the co- 
ordinate point z. The last term in Eq. (6.8) describes the 



change of the density of electrons due to electron-electron 
interactions. 



Note, that due to the condition (6.5), the contribution 



(6.4) of elastic impurities into the collision integral disap- 



pears in Eq. (6.8) after the integration over the direction 
of p. Indeed, Eq. (6.8) is the charge continuity equation 



for the carriers with energy in the interval de, and the 
elastic scatterers cannot affect the corresponding charge 
density. The two terms in the right-hand side of Eq. (6.8 ) 



represent the two inelastic processes changing the density 
{dn^ /de)de: external radiation and electron-electron in- 
teraction. 

The total current through the junction is given by the 
integral of the current density over the energy: 



w 



£F—hQ 



dj]{e,z) dji{e,z) 



de 



de 



de. (6.11) 



Depending on the coordinate point z, it can be carried 
either by the particles in the conduction band or in the 
valence band or both. 

When the junction is non-ballistic, the photocurrent 
depends on the resistances of different parts of the junc- 
tion and of the leads. Let us calculate the photocurrent, 
assuming that the impurities are present only in the "res- 
onant region", \U{z) — erl < while outside this 
region they are absent, and the transport is purely bal- 
listic. The parts of the junction outside this region in 
the leads can be considered as some external circuit, the 
resistance of which can be easily taken into account after 
we obtain the final result for the photocurrent. 



Integrating Eq. (6.8) over the longitudinal coordinate 
z, we obtain 



Sj>(£,zi+0) dj^{e,zi-0) _ 1 



de 



dj^{e,z) 
de 



f ( 



de 
d dn^ 
dt~di 



CewQ, (6.12) 



dz + consti, 



z > zi, (6.13) 



dj^{e,z) 
de 



d dn^ \ 

^ I { I dz^const2, 



z < zi, (6.14) 



where zi{e) is the resonant point, corresponding to the 
energy e. 

Electron-electron collisions conserve the total density 
of particles n{z) = J{dn^/de + dnjde)de at a given 
point z, which leads to the relation 



/ 



de 



d dn^ 

m~di 



d dni 
dt~d^ 



0. (6.15) 



If the ballistic resistance of the junction is much larger 
than the characteristic diffusive resistance induced by im- 
purities and electron-electron interaction in the resonant 
region. 



Rball ^ Ri 



diff^ 



(6.16) 



one can show, that the transmission of particles with en- 
ergy e through the point zo{e), such that U{zo) = e, is 
still strongly impeded, and the photocurrent (6.11 ) in the 



conduction and the valence band is given mainly by the 
currents 



j^is,zi 0) ^^ djiizi+o) shown in Fig.[ll| 



de de ^ 

Indeed, if dn^ is the density of excited electrons m 
the energy interval de in the resonant region, then, for 
instance, the current density 



dj^ {e,zi — 0) dn^ / de 



de 



is much larger than 

dj^{e,zi +0) 
de 



eWRdiffiyo 



dn^/de 



(6.17) 



(6.18) 



since the condition ( 6.16[ ) is fulfilled. 
Taking into account the condition 



de 



_0) 5j>(s,zi+0) 



de 



(6.19) 



and Eq. (6.15), one can find the photocurrent by inte- 



grating the current density at some point z on the left 
from the resonant point: 



w 



dj]{e,z) 
de 

_ ewn^ 

27rv 



de 



(6.20) 
(6.21) 



the same result as given by the Eq. (5.12). 



Thus, we have shown that in a disordered sample, 
where the ballistic considerations of the previous sections 
cannot be immediately applied due to the presence of 
elastic impurities and electron-electron interaction in the 
resonant region, the photocurrent does not change until 
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the disorder or the interactions become too strong, so 
that the impurities-induced or the interactions-induced 
resistance becomes larger than the bahistic one. 

In the opposite hmit, Rdiff ^ Rbaih the excited elec- 
trons and holes, created in the resonant region by the 
radiation, diffuse almost independently of the external 
potential at the p-n interface. The effect of large diffu- 
sive resistance of the resonant region on the photocurrent 
is analogous to the effect of a large resistance of the ex- 
ternal circuit. The photocurrent in this regime is reduced 



rdiff 



Rhall 



diff 



(6.22) 



with Iph given by Eq. (6.21). 

According to the Ref. 46 , in the recent experiments^^"^^' 
the ballistic resistance is of the same order of magnitude 
as the diffusive one, defined as a difference between the 
resistances of n-n and p-n junctions. If Rhaii ^ ^diffi 
Eq. (6.21) is only an order-of-magnitude estimate of the 



photocurrent. As follows from the derivation of the pho- 
tocurrent, the increase of the resistance Rdiff leads to 
the decrease of the photocurrent. One can make the re- 
sistance Rdiff smaller than the ballistic resistance, de- 
creasing the frequency of the EF and thus reducing the 
size of the resonant region and Rdiff- 



VII. POSSIBILITY OF EXPERIMENTAL 
OBSERVATION OF PHOTOCURRENT 

In this section we address the question of experimental 
observability of the photocurrent in a graphene p-n junc- 
tion. We analyze the necessary conditions and estimate 
the value of the photocurrent for achievable radiation in- 
tensities and the junction parameters. 

Geometrical parameters and gate voltages. As dis- 
cussed in the previous section, for the largest photocur- 
rent one needs the diffusive resistance of the junction to 
be smaller than the ballistic one. This can be achieved 
by using sufficiently short junctions^^. Let us take the 
length of the junction L = lOOnm, close to that in 
the experiment in Ref. |12] where the resistance of the 
junction is described by the ballistic model rather than 
by the diffusive one. The typical width W of a p-n 
junction^^^^^^l^^^is a few micrometers; for our estimates 
we take W = Ifim. 

Let Uo = OAeV and Ef = Uo/2 be, respectively, the 
height of the potential barrier and the Fermi energy, close 
to the typical experimental parameters. For the slope of 
the potential we use the naive estimate F = Uq/L, re- 
garding it as constant along the junction. To be more 
precise, considering the effective potential, one should 
take into account the non-uniform charge density distri- 
bution in the junction^"^. However, the corrected in such 
a way potential profile would have the slope of the same 
order of magnitude as the naive estimate. 

Characteristic relaxation lengths. 



For the chosen parameters of a junction the char- 
acteristic length of Landau-Zener tunneling, given by 
Eq. (l3^, is 



A rLz 



13nm. 



(7.1) 



According to the experiment in Ref. where the relax- 
ation time of carriers excited by the near-infrared light 
has been measured, electron-electron interaction is re- 
sponsible for the fastest stage of relaxation, occurring 
on the typical time ^ O.lps, corresponding to the elec- 
tron path lee ^ 0.1/im. Since lee tlz, the condition 
(3.23) is fulfilled, as we assumed in the previous section. 



Another slower stage of relaxation due to the electron- 
phonon interaction has a characteristic length ~ 1/im. 

The mean free path of carriers in graphene is of order of 
1/im at room temperatur^IEl. The characteristic length 
of relaxation due to electron-phonon interaction should 
be of the same order or larger. Then the neglect of such 
a relaxation in the previous section is quite a reasonable 
approximation. 

Desirable radiation frequency. Calculating the current 
due to "2D-modes" in the Sec. [V| we dealt with momen- 
tum scales much larger than = {HF/ttv)^. Accord- 
ingly, to have resonant points on the electron trajectories 
in the "2D-modes" , one should apply the radiation with 
angular frequency ft much larger than 



h 



vF 



(7.2) 



If 1] < Qir^ then the "2D-modes" and some "ID-modes" 
do not have resonant points, that is, there exist electrons 
only weakly affected by the EF and freely penetrating 
through the p-n interface without reflection. As a result, 
if the frequency ft is too low, the photocurrent is strongly 
reduced, and the suppression of tunneling is impossible. 

Note, that according to the Eqs. (5.16) and (5.13) the 
condition Q ^ is equivalent to N2D ^ A^iD, ensuring 
that the junction is effectively two-dimensional. Provided 
this condition is fulfilled, the current is carried mainly by 
the "2D-modes". For our choice of junction parameters 
hQir ^ 29meF, which corresponds to the frequency fir = 
nir/{27r) ^ ITHz. 

As we noted in the Sec. [V| the photocurrent is possi- 
ble if Uq > hVt/2, i.e. hVt < SOOmey (or Vt < 200THz). 
However, to maximize th e ph otocurrent one should sat- 
isfy conditions ( |5.3[ ) and (5.4), for the case under consid- 
eration equivalent to HQ < Uo/2 = 200meV. 

Magnitude of the photocurrent. The characteristic ra- 
diation intensity used in the experiments with nanotube 
junctions^^ is about a few kW/cm'^. Assume, the same 
intensities can be applied to graphene junction, and set 
S = lOkW/cm^^ close to the maximal value reached in 
Ref.ESl 

Then the photocurrent is 



^ ~ 2hcF * 



(7.3) 
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independently of the frequency in the desirable range 
h^lir < hVt < Uo/2. Note, that the photocurrent is a 
few orders of magnitude larger than t hose ob tained in 
the experiments with carbon nanotubeP^'^^'^. 

Possibility of tun nelin g suppression. To maximize the 
dynamical gap, Eq. (3.5), one should lower the frequency. 
For hft = h^lir = 29meV we obtain the dynamical gap 
A ^ 6meV and the exponent of tunneling through it- 
C ^ 10~^, which is insufficient to suppress the tunneling. 
To confine electrons, i.e. to achieve £ 1, one should 
use proportionally larger radiation powers or longer junc- 
tions. 



VIII. CONCLUSION 



To sum up, we studied electron transport in graphene 
junctions irradiated by monochromatic electromagnetic 
field (EF). 

The radiation opens dynamical gaps in the quasiparti- 
cle spectra, proportional to the amplitude of t he E F and 
inversely proportional to its frequency, Eq. (3.5). The 
appearance of the gaps results in a strong modification 
of current- volt age characteristics of a junction. 

If the height of the potential barrier is large enough, 
the directed current {photocurrent), Eq. (5.9), flows 



through the junction without any dc bias voltage ap- 
plied. At small radiation intensities, the photocurrent, 
proportional to the radiation power, Eq. (5.12), is a re- 



sult of inelastic quasiparticle tunneling assisted by one- 
photon absorption. At large intensities, the photocur- 
rent, Eq. (5.15), decreases with radiation power and fi- 



nally saturates at some constant value, Eq. (5.21). 



When the potential barrier is smaller than the pho- 
ton energy h^} but larger than the saturation does 
not happen and the photocurrent decreases to zero at 
large radiation intensities. When the potential barrier 
is smaller than any photocurrent is absent. In 
these regimes, one can adjust the Fermi level in such a 
way, that the quasiparticle transmission in the junction 
is determined by the tunneling through the gap and can 
be fully suppressed, provided that the radiation power is 
large enough. 

In the present paper we also analyze the influence 
of elastic impurities and electron-electron interaction on 
the magnitude of the photocurrent, and show that they 
weakly affect the photocurrent, if the diffusive resistance 
of the junction is not too large compared to the ballistic 
one. 
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APPENDIX A: KINETIC EQUATION IN 
IRRADIATED GRAPHENE 

Now we derive explicitly the kinetic equation, gov- 
erning the dynamics of the electron distribution func- 
tions in graphene exposed to monochromatic electro- 
magnetic wave, taking into account the effect of disor- 
der, electron-electron and electron-phonon interactions 
on the transport. Since graphene in the vicinity of some 
resonant point can be considered as a semiconductor 
with the spectrum, linearized close to the resonant mo- 
mentum, the dynamics of carriers, for which the radi- 
ation matters, should be the same as for conventional 
semiconductor j ^^ l ^^ l ^^ i 

As the radiation can only flip the pseudospin and does 
not induce intervalley scattering, deriving the kinetic 
equation in the lowest non-vanishing order in the radi- 
ation power, we can limit ourselves to the consideration 
of dynamics in a single valley and a single spin direction, 
because the intervalley- and spin- scattering would enter 
only the part of the collision integral, which is indepen- 
dent of the radiation. 

It is convenient to perform calculations in the basis 
of electron states | jp), pseudospin is directed along the 
momentum p, and | |p), pseudospin is antiparallel to 
p. We choose correspondingly the coordinate system in 
the momentum space such that z axis is directed along 
p and the x axis is perpendicular to p and parallel to 
the graphene plane (the plane, in which p can vary), and 
the y axis- normally to the plane. This frame fixes the 
pseudospin basis. 

Since the basis depends explicitly on p, in the momen- 
tum representation one should substitute in the Hamil- 
tonian (3.2) the operator of spatial coordinate r by the 



covariant momentum derivative: 

_ . d Ida 
ap 2 5p 



(Al) 



where a(p) is the angle of rotation of the frame about 
the y axis, normal to the graphene plane. In this section 
h = 1. The second term in the last expression is the 
gauge potential due to the local frame rotations in the 
momentum space. 

The modulus of this term for a given momentum p 
is of order of the corresponding Fermi wavelength Xp ~ 
{da /dp). Far from the Dirac point it is much smaller 
than the characteristic scale on which the potential U (r) 
changes. This allows us to expand the potential up to the 
first order in the gauge field, and write the Hamiltonian 
as 



n ^ vpa, + U{t) + {2p)-^d^U{T)ay 



(A2) 



in the basis chosen. Here dxU{T,p) — {dU{T)/dT)ex{p), 
is the unit vector directed along the perpendicular to 
momentum p the x axis in the graphene plane. 

Now let us proceed to the derivation of the kinetic 
equation for the distribution functions in the basis of 
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states I tp), I Ip). Analogously to the field operator 
4f = (^e,^/^)"^ in a conventional semiconductoi'^, we 
introduce the operator ^ = (^|,^|)-^ with two compo- 
nents acting correspondingly in the conduction and the 
valence bands of graphene. In the momentum represen- 
tation the indices | and | of the operator ^(p) refer to 
the particles with pseudospins aligned along or opposite 
to the momentum p, respectively. Then we introduce the 
non-equilibrium Green's functions 



G>(l,2) = -i(*„(l),*t(2)), 



(A3) 
(A4) 



where a,b =t, i; 1 = {ti,ri}; 2 = {t2,r2}- Accordingly, 
we define^^ the 2x2 matrix Green's functions G"^, G^, 
and G^ , and the matrix function 



G = 



G^ 

G^ 



which satisfies the equatiorP^ 

[(Go ' - S) ® G] = 



(A5) 



(A6) 



(Dyson equation minus its conjugate), where 
GoHl,lO = {idt, - e{l))S{l - V), e is the Hamil- 
tonian of the particles unperturbed by the radiation and 
impurities, and S is the self-energy. The square brackets 
here stand for the commutator [A (g) 5] = B — B ^ A. 
Let us decompose S into two parts. 



s(i,i') = nti)<5(i-i') + s,(i,i'), 



(A7) 



where V{t) is the EF-induced perturbation of the single- 
particle Hamiltonian, 1')- the rest of the self-energy 
part. Assuming that the EF is weak and purposing to 
find the dynamics of the carries in the lowest order in 
the radiation power, we will neglect the effect of external 
radiation on S^. 

It is convenient to solve the problem in the Wigner rep- 
resentation, introducing the "center of mass" coordinates 
T = (^1+^2)72, R = (ri + r2)/2 and the relative ones, 
r = tl — ^2, r = ri — r2, and making Fourier transform 
of all the Green's functions with respect to r and r. 

However, unlike the usual situatioiP^, we expect that 
in the Wigner representation only the quantities 



= 6- vpa, - /7(R) - {2p)-^d,U{Il, p)ay, 



G^^, and are the slow functions of time T, i.e. vary 
on the time scales much larger than the inversed rele- 
vant kinetic energies of electrons, while the other func- 
tions, G^/^ and G^, contain contributions proportional 
to the fast in T perturbation V{T) = -{ev/c)A{T)&. 
The diagonal elements G^ of the Keldysh Green's func- 
tion should be slow since they depend only on the elec- 
tron distribution function, which should vary slowly due 



to the weakness of perturbation. Below we confirm this 
by the direct calculation. 

In order to derive the kinetic equation, we consider the 
Keldysh component of Eq. ( A6 ) : 



[Go^^G^] = [V^G^] 



R ^ riK qK ^ s^A I s^K ^ r^A 



G^-G^0Sf.(A9) 



Taking into account the slowness of G^, one can rewrite 
the left part _of the last equation using the gradient 
approximatio: 



i\ a^G^ -^[/SpG^^ |E„^^,aRG^ 

-i[«pa„G^] + (2p-i)axC/ [<Ty,G^] 



V ®G 



K 



+ T,f ®G^ -G^ ® 
+Ef ®G^-G^®Sf. 



(AlO) 



The brackets [ , ] stand here for the commutators of ma- 
trices, { , }- for the anticommutators. The left part of the 
last equation describes the ballistic properties of elec- 
trons in graphene. The last four terms describe the 
change of distribution functions due to the electron scat- 
tering, independent of the radiation. Further we do not 
consider these terms in detail and focus on the radiation- 
induced transitions, i.e. on the first term in the right- 



hand side of Eq. (AlO). We will only assume, that the 



rate of the radiation-independent scattering between the 
states I tp) and | |p) is not too strong compared to the 
rate of the radiation-induced transitions, so that the lat- 
ter can be calculated independently. The applicability of 
this assumption will be discussed below. 

T he ter m [{dxU /2p)cry^G^] in the left-hand side of 
Eq. (AlO) can be rewritten into the right-hand side (rhs) 
as —{dxU /2p) [o-y (g) G^] in the leading order in the gra- 
dient approximation. Then the term —{dxU /2p)ay can 
be considered in the further calculations as an addi- 
tional small perturbation, induced by the electromag- 
netic fiel d wi th frequency and amplitude going to zero 
[cf. Eq. (3.3)]. Such a perturbation could induce tran- 
sitions between electron states close to the Dirac point. 
However, here we are interested in the radiation-induced 
transitions between the states with large momenta close 
to the resonant ones. Then we will disregard the last 
term in the left-hand side of Eq. (AlO). Within such an 



(A8) approximation the kinetic equation in the form (AlO) in 



the ballistic graphene coincides with that in a conven- 
tional semiconductor, as it should be at large momenta. 

Let us find the EF-induced modification of G^ = G^ + 
G^, disregarding the other processes, the contribution 
of which into the relaxation of the electron distribution, 
as we assumed, can be found separately. Due to the 
radiation some non-zero off-diagonal terms G^ appear. 
In the first order in perturbation in the momentum-time 
representation we obtain 
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G<,{ta,T)Vab{T)G>{T,h)dT 



+ I G>,{ta,T)Vab{T)G<{T,U)dT 



■f 

Jo 



+ / G<,{t,,T)Vat{T)G<{T,h)dT. 

Jta 



(All) 



From now on we assume that in the vicinity of the 
resonance under consideration the RWA can be apphed, 
and that the off-diagonal elements of the perturbation 
are taken in the form 



Vab{T) = Wabe''''^''^ 



(A12) 



Then, from the Eq. ( [All ) we find in the Wigner rep- 
resentation 



G<(T, R,p,e) 



X [ fbS [ e - Eb 



2 



i{ea -Sb^ Qab - iO) 



(A13) 



where we have introduced the total (kinetic+potential) 
energy £a of a particle with momentum p and the pseu- 
dospin a and the distribution function 



/„(r,R,p) = j _G<(r,R,p,£). 

Analogously, 



(A14) 



G>(T, R,p,£) 



i{£a - £6 + ^ab - ^0) 
_ ^ab 

2 



X [{I- fa)5 [e-Sa 
-(l-/5)^fs-£6 + 



(A15) 



Now the obtained functions and ca n be u sed to 
calculate the term [V^G^] in the rhs of Eq. ( |A1Q[ ). The 
convolution of two functions in the Wigner representation 
there can be obtained following the rules formulated in 
the Ref. 53, For instance. 



Vba®G^^ = e-i^^^?VbaGab 



K 



Vba{T)Gi 



^ba 



(A16) 



Before substituting G^ in Eq. (AlO) we integrate over 
the energy e the diagonal element of this equation, in 
order to arrive at the kinetic equation. In that connec- 
tion we also note, that the diagonal elements of the per- 
turbation, Vbb{T) = ^QWbb{^)e'^^^ , do not change the 
distribution functions fa and thus do not enter the ki- 
netic equation. This happens because the longitudinal 



perturbation cannot change the energy level occupation. 
Indeed, 



J de [V^G^]^ 
Y.Wbbme^''^ f dsG^ 



f,(£-^,T) =0. (A17) 



Keeping the second-order terms in EF or the first-order 
ones in the gradients of slow variables, we finally get the 
kinetic equation, 

= 27r\Wba\^S{ea - 65 + ^abXfa ' fb) + (St/)^. (A18) 

Here (St/) 5 is the usual collision integral accounting for 
the change of the distribution function irrespective 
of radiation due to impurity scattering, electron-phonon 
and electron-electron interactions. {crz)bb = (^pI^^^I^p) 
is the direction of the particle velocity in the state with 
momentum p and pseudospin b: (along p or opposite to 

p)- 

If we take the EF-induced perturbation in the form 
of Eq. ( |3.7[ ) within the RWA in the vicinity of the cor- 
responding resonance the kinetics of, e.g., the function 
/l (T, R, p) is described by the equation 



r(p)(/|-/T) + (st/)5, 



where 



r(p) 



27r 



sin^(p,E)(5(2vp-ni]), 



(A19) 



(A20) 



where A is the dynamical gap, and the Planck constant 
is recovered. 



Equation ( A19) can be understood as follows. It differs 



from the usual kinetic equation by the first term in the 
rhs describing the rate of the EF-induced change of the 
distribution function f^. The change is due to the pseu- 
dospin flip under irradiation: | T) ^ I i) and | i) ^ | T)- 
The rate r(p) of both processes can be obtained from the 
Fermi's golden rule. Note, that neither the distribution 
function nor in Eq. ( |A19 ) is assumed to be small; 
the Pauli exclusion principle is taken into account in the 
adduced derivation or in the Fermi's golden rule, so the 
equation is valid for arbitrary distribution functions. 

Let us examine now the applicability of our kinetic 
equation. Deriving the term accounting for the radiation- 
induced transitions, the first term in the rhs of Eq. ( A18), 



we assumed that the distribution functions and /| are 
weakly perturbed by the radiation. As it follows from 
Eqs. (A19) and (A20), this condition is satisfied when 



hvF 



1. 



(A21) 
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Indeed, if we consider the normal incidence of particles 
on a smooth potential barrier, using the obtained kinetic 
equation and setting the collision integral equal to zero, 
we arrive at the tunneling probability 



'- kinetic 



1-A 



(A22) 



which agrees with the result of Eq. (3.19) in the limit 

Another assumption we used when deriving the ki- 
netic equation is that the radiation-independent relax- 
ation processes, such as impurity scattering, electron- 
electron, and electron-phonon interactions weakly influ- 
ence the rate of the radiation-induced transitions that 
can be calculated independently. If we took into account 
this influence, the delta- function in Eq. (A18) would have 



to be substituted by some function smooth on the scale 
^/^Ri where tr is some characteristic time of relaxation 
due to the radiation-independent processes. 

Then, according to Eq. ( A2Q| ), weak radiation would af- 
fect the distribution functions /(T, R, p) in the momenta 
interval of the characteristic width 



5p ' 



h 

trv 



(A23) 



around the resonant Pres = hQ/{2v). For our approxima- 
tion to be valid, the distribution functions in this momen- 
tum interval should be weakly changed by the radiation- 
independent processes. As follows from Eq. (A18), the 



characteristic scale, on which momentum relaxes due to 
these processes, is SpR ~ l/(ri^F), where F is the charac- 
teristic slope of the potential d^U (R) in the region under 
consideration. We need the fulfillment of the condition 
^Pr ^ and obtain thus 



(A24) 



Under this condition the radiation-induce d tran sitions, 
described by the first term in the rhs of Eq. ( A18| ), can be 
considered independently of the other processes described 
by the second term there. 



^ K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 
Y. Zhang, S. V. Dubonos, L V. Grigorieva, and A. A. 
Firsov, Science 306, 666 (2004). 

2 P. R. Wallace, Phys. Rev. 71, 622 (1947). 

^ A. K. Geim and K. S. Novoselov, Nature Mat. 6, 183 
(2007). 

^ S. V. Morozov, K. S. Novoselov, F. Schedin, D. Jiang, 
A. A. Firsov, and A. K. Geim, Phys. Rev. B 72, 201401 (R) 
(2005). 

^ K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 
M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, and 
A. A. Firsov, Nature 438, 197 (2005). 

^ Y. Zhang, Y.-W. Tan, H. L. Stormer, and P. Kim, Nature 
438 (2005). 

^ A. Das, S. Pisana, B. Chakraborty, S. Piscanec, S. K. Saha, 
U. V. Waghmare, K. S. Novoselov, H. R. Krishnamurthy, 
A. K. Geim, A. C. Ferrari, et al., Nat. Nanotechnol. 3, 201 
(2007). 

^ S. V. Morozov, K. S. Novoselov, M. I. Katsnelson, 
F. Schedin, D. C. Elias, J. A. Jaszczak, and A. K. Geim, 
Phys. Rev. Lett. 100, 016602 (2008). 

^ Y. Q. Wu, P. D. Ye, M. Capano, Y. Xuan, Y. Sui, M. Qi, 
and J. Cooper, Appl. Phys. Lett. 92, 092192 (2008). 
K. L Bolotin, K. J. Sikes, Z. Jiang, G. Fudenberga, 
J. Hone, P. Kim, and H. L. Stormer, arXiv e-prints (2008), 
0802.2389. 

M. C. Lemme, T. J. Echtermeyer, M. Bans, and H. Kurz, 
IEEE Electron Device Lett. 28, 282 (2007). 
^2 B. Huard, J. A. Sulpizio, N. Stander, K. Todd, B. Yang, 
and D. Goldhaber-Gordon, Phys. Rev. Lett. 98, 236803 
(2007). 

J. R. Williams, L. DiCarlo, and C. M. Marcus, Science 
317, 638 (2007). 



B. Ozyilmaz, P. Jarillo-Herrero, D. Efetov, D. A. Abanin, 
L. S. Levitov, and P. Kim, Phys. Rev. Lett. 99, 166804 
(2007). 

Z. Chen, Y.-M. Lin, M. J. Rooks, and P. Avouris, Physica 
E 40, 228 (2007). 

R. V. Gorbachev, A. S. Morozov, A. K. Savchenko, D. W. 
Horsell, and F. Guinea, arXiv e-prints (2008), 0804.2081. 
G. Liu, J. Valesco Jr., W. Bao, and C. N. Lau, arXiv e- 
prints (2008), 0804.2513. 

G. W. Semenoff, Phys. Rev. Lett. 53, 2449 (1984). 
F. D. M. Haldane, Phys. Rev. Lett. 61, 2015 (1988). 

M. L Katsnelson, K. S. Novoselov, and A. K. Geim, Nature 
Phys. 2, 620 (2006). 

X. Duan, Y. Huang, Y. Cui, J. Wang, and C. M. Lieber, 
Nature 409, 66 (2001). 

M. S. Marcus, J. M. Simmons, O. M. Castellini, R. J. Ham- 
mers, and M. A. Eriksson, J. Appl. Phys. 100, 084306 
(2006). 

M. Freitag, Y. Martin, J. A. Misewich, R. Martel, and 

P. Avouris, Nano Lett. 3, 1067 (2003). 

Y. Ohno, S. Kishimoto, T. Mizutani, T. Okazaki, and 

H. Shinohara, Appl. Phys. Lett. 84, 1368 (2004). 

J. Wang, M. S. Gudiksen, X. Duan, Y. Cui, and C. M. 
Lieber, Science 293, 1455 (2001). 

P. Hanggi, in Quantum transport and dissipation, edited by 
T. Dittrich, P. Hanggi, G.-L. Ingold, B. Kramer, G. Schon, 
and W. Zwerger (Wiley- VCH, Weinheim, 1998). 

^'^ V. M. Galitskii, S. P. Goreslavskii, and V. F. Elesin, Sov. 
Phys. JETP 30, 117 (1970). 

2^ V. F. Elesin, Sov. Phys. Semicond. 4, 1302 (1970). 

2^ V. F. Elesin, Sov. Phys. Semicond. 11, 1470 (1970). 

2° V. F. Elesin, Sov. Phys. JETP 32, 328 (1971). 

V. F. Elesin, A. L Erko, and A. L Larkin, JETP Lett. 29, 



21 



651 (1979). 

^2 Q. T. Vu, H. Haug, O. D. Miicke, T. Tritschler, M. We- 
gener, G. Khitrova, and H. M. Gibbs, Phys. Rev. Lett. 92, 
217403 (2004). 

S. P. Goreslavskii and V. F. Elesin, JETP Lett. 10, 316 
(1969). 

A. S. Alexandrov and V. F. Elesin, Sov. Phys. JETP 35, 
403 (1972). 

T. J. Inagaki and M. Aihara, Phys. Rev. B 66, 075204 
(2002), and references therein. 

M. V. Fistul and K. B. Efetov, Phys. Rev. Lett. 98, 256803 

(2007) . 

M. Moskalets and M. Biittiker, Phys. Rev. B 66, 205320 
(2002). 

E. J. Mele, P. Krai, and D. Tomanek, Phys. Rev. B 61, 
7669 (2000). 

L. Allen and J. H. Eberly, Optical resonance and two-level 
atoms (Dover, New- York, 1975). 

E. O. Kane and E. L Blount, in Tunneling Phenomena in 
Solids, edited by E. Burstein and S. Lundqvist (Plenum 

Press, New York, 1969). 

V. V. Cheianov and V. L Fal'ko, Phys. Rev. B 74, 
041403(R) (2006). 

A. V. Shytov, N. Gu, and L. S. Levitov, arXiv e-prints 

(2008) , 0708.3081. 

M. V. Fistul and K. B. Efetov, Phys. Rev. B 76, 195329 

(2007) . 

L. M. Zhang and M. M. Fogler, Phys. Rev. Lett. 100, 
116804 (2008). 

M. Abramowitz and L A. Stegun, eds.. Handbook of Mathe- 
matical Functions with Formulas, Graphs, and Mathemat- 
ical Tables (U.S. Government Printing Office, 1972), 10th 
ed. 

M. M. Fogler, D. S. Novikov, L. L Glazman, and B. L 
Shklovskii, Phys. Rev. B 77, 075420 (2008). 
^'^ L L. Aleiner and K. B. Efetov, Phys. Rev. Lett. 97, 236801 
(2006). 

E. G. Mishchenko, Phys. Rev. Lett. 98, 216801 (2007). 

J. M. Dawlaty, S. Shivaraman, M. Chandrashekhar, 

F. Rana, and M. G. Spencer, Appl. Phys. Lett. 92, 042116 

(2008) . 

^° K. Henneberger, Physica A 150, 419 (1988). 

K. Henneberger and H. Haug, Phys. Rev. B 38, 9759 
(1988). 

F. Jahnke and K. Henneberger, Phys. Rev. B 45, 4077 
(1992). 

J. Rammer and H. Smith, Rev. Mod. Phys. 58, 323 (1986). 



